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Based on Sjogren und Sjolander's Mode-Coupling(MC)-Model, we have reformulated and calculated 
the memory kernel (MK) of the velocity autocorrelation function (VACF)on a Nio.2Zro.g-liquid. 
Reformulating means here that we have constructed the memory kernel of VACF for our binary 
system, instead of one for one atomic system of the Sjogren und Sjolander's model. The data required 
for the theoretical calculations have been obtained from molecular dynamics (MD) simulations. The 
theoretical results then are compared with those directly obtained from computer simulation. We 
found, although it exists a qualitative agreement between theoretical predictions and simulation 
results, that quantitatively there is an deviations between both results, especially for Zr-subsystem. 
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I. INTRODUCTION 

^— > 

Since last decade and until today, one important goal of modern research in condensed matter is to find a good 
quantitative description of the glassy dynamics in liquids. Theoretically there is a established theory, the so-called 
mode-coupling theory (MCT), that can describe and predict important parameters and quantities in glassy systems 
qualitatively good. Work in the last decade has provided evidence pj that MCT is able to describe the slow 
dynamics of fragile liquids in the weakly supercooled state. The recent work of Kob et.al.|3| has also showed that 
theory able to accurately describe the non-ergodicity parameters of simple as well as of network- forming liquids. The 
one important statement of the MCT is that the dynamics processes at long times is driven by a memory kernel which 
is included in MCT's integro-differential equation. 

In previously paper Q we have studied the memory kernel of incoherent intermediate scattering function. Further, 
we have compared the memory kernels Mf(q,t) evaluated from the MCT formula with the kernels M®(q,t) from 
, inverting the time evolution of the intermediate scattering functions. The comparison shows encouraging agreement 
\Q ' at 900 to 1100 K while significant deviations are found at 1200 K. It is an open question whether for this temperature 
an improved agreement between calculated Mf(q,t) and estimated Mf(q,t) can be obtained by inclusion of the 
coupling to transversal currents in the memory kernel formula as provided, e.g., by Gudowski et.al. 

The main objective of paper is to give one of the answer of above open question. But instead of the memory kernel 
of incoherent intermediate scattering function, we investigated the memory kernel K(t) of the velocity autocorrelation 
function (VACF) for temperature T = 1500 K. By calculating the memory kernel we included all coupling as proposed 
by Sjogren und Sjolander's Mode-Coupling(MC)-Model The calculated memory function then is compared with 
the kernel from inverting the time evolution of VACF. 

Our paper is organized as follows: In Section lllll we present the model and give some details of the computations. 
Section ITT1 gives a brief discussion of some aspects of the Sjogren und Sjolander's Mode-Coupling(MC)-Model as used 
£h ! here. Results of our MD-simulations and their analysis are then presented and discussed in Section llVl 
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II. SJOGREN-SJOLANDER'S MODE COUPLING MODEL 

As we know, the VACF can be studied through the formalism developed by Zwanzig and Mori [T^j > which is based 
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C3 on a following integral equation: 

dm = r* 

dt io 



K{t')ip{t - r)dt' (1) 



In this paper the Sjogren und Sjolander's (SS) Mode-Coupling(MC)-Model [|| is used. In this model memory 
kernel K (t) of VACF ip(t) was proposed for a simple liquid according to a combination of kinetic and mode coupling 
concepts. 

The basic idea of the model is that memory-kernel can be divided into two term. The first term comes from the 
uncorrelated binray collisions and the other from correlated collisions. The former is related to the fast decay of kernel 
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at short times. The latter is expresses the mode-coupling term which incorporate more sophisticated processes that 
appear at longer times. Moreover, this term is based on the idea that the motion of a tagged particle is influenced 
by constraints collectively imposed by its neighboring particles. The memory kernel K(t) of VACF take also the 
following relation 



K(t)=K B (t)+K MC (t) 



(2) 



According to SS-Model the mode-coupling effects take into account four different coupling-term: density-density 
coupling, two contributions from density-longitudinal current couplings, and density-transverse current coupling (s. 
Sjogren (1980/81) dEl, Gudowski et.al. (1993) 0, Canales et.al. (1997) [13). 

To compute the memory kernel of VACF, we follow method used by Sjogren (1980/81), Gudowski et.al. (1993), 
Canales et.al. (1997). If by these authors one atomic systems only were investigated, here we investigate a binary 
system. Also we have to develop and construct the formulation of coupling terms for our binary system. 

Here are the results of our reformulation of K(t) for our binary system. The binary term is defined by |l4j 



where £l% a is the Einstein-frequency of a-component ^ 
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and T a can be determined by the following equation 
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K(0) 



(3) 



(4) 



(5) 



The coupling terms is expressed by 

Kl IK (t) = K 0Qa (t) + K 01a (t) + K lla (t) + K 22a {t). 



(6) 



Here Koo a (t) describes the density-density coupling, Koi a (t) and Kn a (t) the density-longitudinal currents couplings, 
and K 22a (t) the density-transverse current coupling. 
The Laplace transforms of these terms can be written 



(7) 



K 01a (z) = Roi a (z) K Ba (z) + K a {z) 



Kn a (z) = K Ba (z)Rii a (z)K a (z), 



(8) 
(9) 



K 2 2a = K Ba (z) + Raoa(z) + K Ba (z)R Qla (z) R 2 2 a (z)K a (z) 



(10) 



where K Ba (z) und K a {z) are repectively the Laplace transform of the binary term and the total memory kernel. The 
Laplace transform of the total Memory kernel is as follows 



K a (z) = K Ba (z) + R QOa (z) + K Ba (z)R 01a (z) /X 



(11) 



with 



X = 1 - Rai a (z) - K Ba {z)Rx\ a {z) - 

K Ba (z) + R 00a (z) + K Ba (z)Roi a (z) R 22a (z) 



(12) 



3 



where Rij a are the Laplace transform of " recollision" terms Rij a (t). Assuming isotropic systems Rij a (t) can be 
expressed as follows 



with 
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q 4 (cF) a (q,t)AF sa (q,t)dq, 



(13) 



(cF) a (q,t) = [c aa {q)f F aa (q,t) + 

2c aa (q)c Q p(q)F a 0(q,t) + [c a(3 (q)} 2 F aa (q,t); a ^ 



(14) 



with 



Roia(t) 



q 2 (LF') a (q,t)AF sa {q,t)dq, 



(15) 



(LF r ) a (q,t) = c aa (q) 
Caa(q) 



k B Tq 2 
IhaaW H c aa {q) 

777 X qi 

k B Tq 2 ' 

iLapyq) H c a p(q) 

m a Xf} 

I \ . k B Tq 2 ' 

jLaaW H C Qct (g) 

777 q, iZ- q< 

iLapW H C Q/3 ((j) 



(</,*) + 
^ (<?,*) + 
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with 



oo 

fliiaW = ~ fi 2 o4 / q 2 {BC L ) a (q,t)AF sa (q,t)dq, 



(BC L ) a (q,t) 



I , . k B Tq 2 



777 Q' X fy^ 



I s , k B Tq 2 
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C Laa {q,t) + 



k B Tq 2 
Ihapyq) H c Q/3 (g) 



k B Tq 

lLap{q) + c a p{q) 

m a x (i 
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(17) 



(18) 



with 



oo 

7W*) = - R 2 i / q 2 (TC T ) a (q,t)AF sa (q,t)dq, 



(19) 



(TC T ) a (q,t) = [lTaa{q)f C Toia (q,t) + 

2iT aa (qhT a i3(q)C Ta f3(q,t) + [iT a p{q)f C Taa {q,t) 



(20) 



Here c Qa (g) is the Fourier transform of aa— direct correlation function, ^Laa(q) and 7tcm(<z) are q -dependent 
quantities, which are defined in Balucani |14| . Cr cta (g,i) und CL aa (q,t) are the transverse and longitudinal current 
correlation functions of the aa- parts. AF sa (q,t) is the difference between F sa (q,t) and F 0a (q,t), where F 0a (q,t) — 
exp(— (k B T /2m a )q 2 t 2 ) expresses the free-particle-form of intermediate incoherent scattering function. 
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III. SIMULATIONS 

The simulations are carried out as state-of-the-art isothermal-isobaric (N,T,p) calculations. The Newtonian equa- 
tions of TV = 648 atoms (130 Ni and 518 Zr) are numerically integrated by a fifth order predictor-corrector algorithm 
with time step At = 2.5 10 _15 s in a cubic volume with periodic boundary conditions and variable box length L. With 
regard to the electron theoretical description of the interatomic potentials in transition metal alloys by Hausleitner 
and Hafner Q, we model the interatomic couplings as in 8] by a volume dependent electron-gas term E vo i(V) and 
pair potentials <p(r) adapted to the equilibrium distance, depth, width, and zero of the Hausleitner-Hafner potentials 
Q for Ni 2 oZr 8 o For this model, simulations were started through heating a starting configuration up to 2000 K 
which leads to a homogeneous liquid state. The system then is cooled continuously to various annealing temperatures 
with cooling rate —dtT = 1.5 10 12 K/s. Afterwards the obtained configurations at various annealing temperatures 
(here 1500-800 K) are relaxed by carrying out additional isothermal annealing runs. Finally the time evolution of 
these relaxed configurations is modelled and analyzed. More details of the simulations are given in [jj . 

IV. RESULT AND DISCUSSION 

We have calculated the total memory kernel K a (t) of VACF in range 0, 024 < q < 4, 796 A -1 . The input data are 
all correlation functions that are analyzed from MD data. For simplifying we have only calculated the memory kernel 
at T = 1500 K. The reason for this is, because (i) according to some authors, except binary term and density-density 
coupling, the other couplings could be neglected by lower temperature, i.e., the contribution of the other couplings 
are very small, and more smaller than the contribution of the binary term and density-density coupling 16], (ii) after 
our experience there is a coupling contribution which decays slowly at lower temperature. These slowly decays make 
us technically a problem by the Laplace transform of couplings. 

The figures an d © shows our results of the memory kernel K a (t) of VACF. MD-results can be computed over 
the equation Q). From figures we can see that the results of SS-model do not have a good agreement with MD-results, 
escpecially those for Zr-atom. By Zr-atom the run of the memory kernel K a (t) shows some of oscillations that by 
MD-results do not appear. These oscilations are primaly from the density-density coupling. According to ea. ll4|l 
the intermediate coherent scattering functions F(q,t) are responsible for those oscilations. The deviations between 
SS-models results and MD-results can also be shown by results in computing tp(t), e.g., on figures 10 and (@J. 

As shown in figures Q and J2J the most significant contributions come from the density- density coupling, the second 
from the density-longitudinal current density coupling, and the least from the density-transverse current coupling. 
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FIG. 1: Ni - K(t)/K{0). 



V. CONCLUDING REMARKS 

We have calculated the memory kernel of VACF based on SS-model. The parameters used by calculating have 
produced from MD-simulation. We have compared the SS-model's results and MD-results. Both results are shown a 
qualitative agreement, although quantitatively there is a deviations between both results, especially for Zr-subsystem. 
These results are agreed with the generall Mode Coupling Theory's predictions 0, namely that the most significant 
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FIG. 2: Zr - K(t)/K(0). 
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FIG. 3: Ni - tP(t). 

contributions are resulted from the density-density coupling. This coupling always is dominant when the temperature 
of system is lowered, and especially when the temperature of system is near the critical glass temperature T c . 
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FIG. 4: Zr - ip(t). 
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